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Abstract 



> . 

We consider a nonholonomic system describing a rolling of a dynamically non- 
' symmetric sphere over a fixed sphere without slipping. The system generalizes the 

. classical nonholonomic Chaplygin sphere problem and it is shown to be integrable for 

• ' one special ratio of radii of the spheres. After a time reparameterization the system 

becomes a Hamiltonian one and admits a separation of variables and reduction to Abel- 
qq ' Jacobi quadratures. The separating variables that we found appear to be a non-trivial 

C^) . generalization of ellipsoidal (spheroconic) coordinates on the Poisson sphere, which can 

be useful in other integrable problems. 

Using the quadratures we also perform an explicit integration of the problem in 
' theta-functions of the new time. 

1 Introduction 

One of the best known integrable systems of the classical nonholonomic mechanics is 
the Chaplygin problem on a non-homogeneous sphere rolling over a horizontal plane 
without slipping. In [Sj S. A. Chaplygin obtained the equations of motion, proved 
their integrability and performed their reduction to quadratures by using spheroconical 
coordinates on the Poisson sphere as separating variables. He also actually solved the 
reconstruction problem by describing the motion of the sphere on the plane. 

Various aspects of this celebrated system were studied in [H[T9j[8][20] > anc ^ ^ s ex P u °it 
integration in terms of theta-functions was presented in [17] . 
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Several nontrivial integrable generalizations of this problem were indicated by V. 
Kozlov [21] (the motion of the sphere in a quadratic potential field), A. Markeev [23] 
(the sphere carries a rotator), in [25| (an extra nonholonomic constraint is added) and in 
[15] (the sphere touches an arbitrary number of symmetric spheres with fixed centers) . 

Next, amongst others, the papers [3J HI considered rolling of the Chaplygin (i.e., 
dynamically non-symmetric) sphere over a fixed sphere, so called sphere-sphere problem. 
They studied the equations of motion in the frame attached to the body. More generic 
(although more tedious) form of the equations also appeared in the works of Woronetz 
[571 [55], who, nevertheless, solved a series of interesting problems describing rolling of 
bodies of revolution or flat bodies over a sphere. 

A rolling of a generic convex body over a sphere was also discussed in the recent 
survey [5]. 

In O HI [26] it was observed that the equations of motion of the Chaplygin sphere- 
sphere problem admit an invariant measure, but, as numerical computations show, in 
the general case they are not integrable (there is no analog of the linear momentum 
integral). However, as was also found in [3], for one special ratio of radii of the two 
spheres an analog of such an integral does exist, and the system in integrable by the 
Jacobi last multiplier theorem. 

Until recently, no one of the above generalizations was integrated in quadratures 
(except a case of particular initial conditions of the Kolzov generalization, considered 
in PI]). 

In this connection it should be noted that a Lax pair with a spectral parameter for 
the Chaplygin sphere problem or for its generalizations is still unknown and, probably 
does not exists. Hence, one cannot use the powerful method of Baker- Akhieser functions 
to find thcta-function solution of the problem. 

Contents of the paper. Our main purpose is to find appropriate separating vari- 
ables, which allow to reduce the integrable case of the Chaplygin sphere-sphere problem 
to quadratures, as well as to give explicit theta-function solution. 

It appears that, in contrast to the classical Chaplygin sphere problem, the usual 
spheroconical coordinates on the Poisson sphere do not provide separating variables 
and that such variables should be introduced in a more complicated way (see formulas 

Using the quadratures, we also give a brief analysis of possible bifurcations and 
periodic solutions. These results are presented in Sections 3-4. 

In Section 5 we briefly describe another type of periodic solutions. 

Section 6 provides a derivation of explicit thcta-function solutions of the problem in 
a self-contained form f Theorem 16. 4p . 

Finally, in Appendix we show how the separating variables we used can be obtained 
in a systematic way, by reducing a restriction of our system to an integrable Hamilto- 
nian system with 2 degrees of freedom and applying a classical result of Eisenhart on 
transformation of the Hamiltonian to a Stackel form. 

2 Equations of motion and first integrals 

Consider rolling of the Chaplygin sphere inside/over a fixed sphere without slipping 
(Figure 1.) 

Let w, m, I = diag(/i, 12, I3), and b denote respectively the angular velocity vector 
of the Chaplygin sphere, the mass of the sphere, its inertia tensor, and the radius. 
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By n we denote the unit normal vector to the fixed sphere S 2 at the contact point 
P. The angular momentum M of the moving sphere with respect to P can be written 

as 

M = Iiv + dn x (n x u), d = mb 2 , (2.1) 

where x denots the vector product in R 3 . 

The phase space of the dynamical system is the tangent bundle T(SO(3) x S 2 ). By 
using the no slip nonholonomic constraint (which corresponds to zero velocity in the 
point of contact), one can obtain the reduced equations of motion in the frame attached 
to the sphere in the following closed form (see, e.g., [3]): 

M = Mxlj, n = knxu>, k = — — -, (2.2) 

a + b y ' 

a being the radius of the fixed sphere. 

Note that the ratio k can take any positive or negative value depending on the 
relative position of the rolling and fixed spheres, as illustrated in Fig. 1. 




n 



a) < k < 1 




b) k > 1 c ) k < 



Figure 2.1: Rolling of the Chaplygin sphere inside/over the fixed sphere (dashed). 

For arbitrary k the equations (|2.2p possess three independent integrals 

Fo = (n,n) = l, H={M,u), F 1 = (M,M), (2.3) 

which, in view of (|2.ip , can be written as 

F = (n,n) = l, H=(u>,Ju J )-d 2 (n,u>) 2 , 
Fi = (Ju>, Jw) — 2(Ju>, n) (n, u>) + d (n, u>) , 

where J = I + rfE, E being the identity matrix. 
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As shown in [26] , the equations (|2.2p expressed in terms of omega, n also have the 
invariant measure p du) dn with the density 

p = y/(n, n) - d{n, J~ l n). (2.5) 

Thus, according to the Jacobi theorem, for a complete integrability of this system one 
extra integral is needed. 



The Chaplygin sphere on the plane. Clearly, the case k = 1 corresponds to 
a — ► oo, that is, the fixed sphere transforms to a horizontal plane with the unit normal 
vector n, and we arrive at the classical integrable Chaplygin problem, when the linear 
(in M) momentum integral is preserved: 

(M,n) = (Iu},n). (2.6) 



Second integrable case. According to [3J, the system (|2.2j) is also integrable in 
the case k — — 1, which describes rolling of a non- homogeneous ball with a spherical 
cavity over a fixed sphere and the quotient of the radii of the spheres equals - = 5 (see 
Fig. 1 c). 

In this case, instead of (|2.6[) . there is the following linear integral 

F 2 = (AM, n) , (2.7) 

where 

A = diag(J 2 + J 3 - Ji, J3 + J1- J2, J1 + J2- J3) ■ 

Note that, as was shown in [BJ, the modification of this system obtained by imposing 
the extra "no twist" constraint (w, n) = (sometimes called as the rubber Chaplygin 
ball) is also integrable for the ratio k = — 1. 

In the next sections we present explicit integration of this case under the condition 
F2 = 0. Our procedure is similar to that of the problem of the Chaplygin sphere rolling 
on a horizontal plane in the case of zero value of the area integral (|2.6p , (see [HI EH [7] ) , 
however, analytically, it is more complicated. 

Integration of the system (|2.2p with k — — 1 in the general case F 2 7^ is still an 
open problem. 



A remark on reduction to quadratures in the case d = 0. Note that in 
the limit case d — one has M = Iiv and the equations (|2.2p with k = — 1 take the 
form 

la; = lu) x ui, n = —n x u>. 

As was noticed in [3], by the substitution 

M = Ala;, 7 = n (2.8) 

and the sign change t — * — t, the latter system transforms to the Euler-Poisson equations 
for the classical Euler top, 

Ai = M. x w, 7 = 7x0;, (2.9) 
uj t . = a t Mi, a l = — — (2.10) 



1 We could not interpret this integral as a momentum conservation law. 
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which possesses first integrals 



(At,7> = fl, (M,AM) = h, (M,M) = f, 

where A. = diag(ai, a 2 , as). 

As was indicated in several publications (see, e.g., [HI]), the Euler-Poisson equations 
(|2.9p can be integrated by separation of variables. Namely, by an appropriate choice 
of the constant vector 7 in space we can always set g = 0. Then (.4a;, 7) = 0, and 
the equations (|2.9p reduce to a flow on the tangent bundle of the Poisson sphere S 2 — 
{(7j7) — !}■ I n the spheroconical coordinates Ai, A2 on S 2 such that 

7!= (a i -X 1 )(a i -X 2 ) (2<u) 
(^ - a,-)(aj - a*,) 



the flow is reduced to the quadratures 



Xi dX\ X2 (IX2 



= 0, 

(2.12) 

= C dt, C — const, 



R(X) = -(A - oi)(A - oa)(A - 03) (/A - h). 

The latter contain one holomorphic and one meromorphic differential on the elliptic 
curve £ = {^i 2 = R(X)}. Thus the the quadratures give rise to a generalized Abel-Jacobi 
map and, following the methods developed in [10j . they can be inverted to express the 
variables -y,u> in terms of theta-functions of £ and exponents (see, e.g., [THUS] for the 
concrete expressions). 

Apparently, in the general case d ^ the substitution (|2.8p is not useful to integrate 
the the equations (|2.2[) in the second integrable case fc = — 1. In particular, it does not 
transform these equations to the case of the classical Chaplygin sphere problem (fc = 1). 



3 Reduction to quadratures in the case F2 = 

We now consider the case d ^ 0, but assume that the linear integral F2 in (|2.7p is zero, 
which imposes restrictions of the initial conditions. Then, from (|2.2p with k = — 1 and 
F2 — we get 

n = n x u>, (u>, Bn) = 0, 

where B = ( J — dn <8> n)A. This allows to express the angular velocity in terms of fi, n 
in the following homogeneous form 

Bn x fi , 
(n, Bn) v ; 

In view of the above remark on the reduction to quadratures in the case d = 0, to 
perform separation of variables it seems natural to use the spheroconical coordinates 
Ai, A2 given by (|2.1ip (with 7, replaced by m) in the general case d 7^ too. However, 
this choice does not lead to success: after some calculations one can see that the first 
integrals H,F% have mixed terms in the derivatives Ai, A2. 

It appears that a correct choice is given by the following quasi- spheroconical coordi- 
nates z\, z 2 on the Poisson sphere (n,n) = 1: 

<= r( l m dC .^ 7 ( M ,fc). (1,2,3), (3.2) 

G(zi, z 2 ) (Ji ~ djJjJk (ai - a J )(a l - a k ) 
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where 



G(z 1 ,z 2 ) = 1 - d(TrJ - 2d)(z 1 + z 2 ) + d(4dct J — dTr(JA))ziz 2 , (3.3) 

and, as in (|2.10p , a, = (AiJi)" 1 . A systematic derivation of the substitution (|3.2|) is 
presented in Appendix. 

Note that when d = 0, the factor G becomes 1 and the relation (|3.2p takes the form 
of (|2.11|) . that is, zi, z 2 do become the usual spheroconical coordinates on S 2 . 

We note that similar quasi-spheroconical coordinates were already used in [6j [5] to 
integrate the "rubber" Chaplygin sphere-sphere problem. 

In the above coordinates z\ , z 2 one has 



p 2 = ((In, J" l n) = fl\ r( l y <n,Bn) = detldetA-^ 
detJG(zi,Z2j G(zi,z 2 ) 



and 



z\ 



Z-2 



G(z 1 , z 2 ) 

2 \ z\ - Oi z 2 -a,i G(z!,z 2 ) 

Then the expressions (|3.1|) yield 



7&j fife J j Jfc 

2 Ji — d 



(1 - dA i z 2 )z 1 



(1 - dAiZi)z 2 



(1-a zi)(l-a fe zi)z 2 (l-o,- z 2 )(l - a fe z 2 )zi 



(u>, n) 



n\n 2 n z (Ji - J 2 )(J 2 - J3XJ3 ~ Ji) z 2 ) 
2 dct / 
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Z-2 



<S>{z 1 )z 2 $(z 2 )zi 



(3.4) 
(3.5) 

(3.6) 
(3.7) 



$(z) = (a^z - l)(a 2 " 1 z - l)(a 3 " 1 z - 1). 



Substituting them, as well as (|3.2p . into the integrals H, F\ in (|2.4|1 . after simplifications 
we get 



H=( Zl - z 2 ) 
Fx = ( Zl - z 2 ) 



dct I 



4G 2 (zi,z 2 ) 
detl 

4G 2 (zi,z 2 ) 



2 *2 



gfo) , 2 £(fl) .; 

9f{zi)4 1 $(z 2 )z 2 

-0(^2) . 2 1p{zi) 



$(zi)» 



2*1 



*(Z 2 )-Z 



2*2 



(3.8) 



where 

*(z) = ddetAz 2 -Tr(AJ)z + 2, V(z) = (4 dot J — dTr(AJ))z — (TrJ — 2d). (3.9) 
Next, substituting (|3~B)) . (j3~7| . and (|3~2)) into (|2~Tj) , we also obtain 

Zl Z2 



(1 - a j 1 z 1 )(l - a k 1 zi)z 2 (1 - a j 1 z 2 )(l - a k 1 z 2 )z 1 



(3.10) 



Now, fixing the values of the integrals by setting H = h, F± = f , then solving 
with respect to if, z\ and using the relation 



we get 



*(z 2 M^i) - y(zi)ip{z 2 ) = det A(z 2 - z 1 )G{z ll z 2 ), 



zl G{z!,z 2 ) 4$(z a )(/*(z a ) + H{z a )) 



(3.H) 



(zi - z 2 ) 2 detldet A 

,) 

2 A(z a - ai)(z a - a 2 ){z a - a 3 ) (f&(z a ) + hip{z a )), (a, (3 = (1, 2). 



G(z!, Z 2 ) 
(Zl - Z 2 ) 2 
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After the time reparameterization 

1 



dt = 



y/G(zi,Z 2 ) 



-.dT 



det J 
dctl 



(In, J 1 n) dr 



the above relations give 



dz\ 
~dr~ 



Z 2 



dz 2 

~d7 



*2 



Zl - Z 2 CLT Z 2 - Zl 

R{z) = -{z - ai)(z - a 2 )(z - a 3 ) (f^(z) + H{z)). 
The latter are equivalent to the following Abel-Jacobi type quadratures 

dz\ dz 2 zi dz\ z 2 dz 2 



2dr, 



0. 



(3.12) 



(3.13) 



(3.14) 



which contain 2 holomorphic differentials on the hyperelliptic genus 2 curve T = {w 2 = 
R(z)}. 

In view of (|3.9p . when d — > 0, the polynomial R(z) becomes a degree 4 polynomial, 
and (|3.14|) reduce to the quadratures (|2.12|) for the Euler top problem, as expected. 

It is interesting that, like in the integration of the original Chaplygin sphere problem 
presented in [8] , the reparameterization factor in (|3.12[) coincides with the density (|2.5[) 
of the invariant measure. In the real motion this factor never vanishes, hence the 
reparameterization is non-singular. 

Now substituting the above formulas for if, if m t° fl3.10|l . ()3.6|1 and simplifying, 
we express the angular momentum M and the velocity u> in terms of z\,z 2 and the 
conjugated coordinates wi,w 2 : 



Mi = 



y/(aj - Zi)(aj - Z 2 ) \/(a k - zi)(a k - z 2 ) 



y/iJi-JjXJi-Jk) 



zi - z 2 



Wl 



2y/G( Zl ,Z 2 ) 

W 2 



(zi - Oj)(zi - afc) (z 2 - a,j)(z2 - ait) 



\fljh 



y/(a,j - z 1 )(a j - z 2 ) y/(a k - z 1 )(a k - z 2 ) 



y/(Ji- Jj)(Ji- Jk) 



Zl ~ Z 2 



2^G(z 1; z 2 ) 

(gi ~ l/(dA z ))w 2 

(zi - Oj)(zi - a k ) (z 2 - cij)(z 2 - a k ) 



(z 2 - l/(dAi)wi 



(3.15) 



(3.16) 



Next, the projection (u),n) in f|3 . 7|) takes the form 

(w,n) = JiJ 2 J 3 V (zi - ai)(zi - a 2 )(zi - a 3 ) (z 2 - ai)(z 2 - a 2 )(z 2 - o 3 ) 



zi - z 2 



Wi 



w 2 



(zi - ai)(zi - a 2 ){zi - a 3 ) (z 2 - ai)(z 2 - a 2 )(z 2 - a 3 ) 



Since 



u a = \J-[z a - ai)(z a - a 2 )(z a - a 3 ) (/f(2 a ) + hip(z a )) 



= \/-{z a - a{)(z a - a 2 )(z a - a 3 ) fddetA(z a - ci)(z a - c 2 ), 
the latter relation also reads 



(u>,ri) = JiJ 2 J 3 y/{a - zi)(a - z 2 ) y/ (c 2 - zi)(c 2 - z 2 ) 



■1 



Zl - Z 2 



Wi 



w 2 



(zi - ci)(zi - c 2 ) (z 2 - ci)(z 2 - c 2 )_ 

In Section [S] we shall use the above expressions to obtain explicit theta-function 
solutions for the components of n, to, M in terms of the new time r. 
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4 Qualitative study of the motion and bifurcations 



For generic constants h, f the polynomial R(z) in (|3.14[) has simple roots dj, c\, c 2 , and 
in the real motion the separating variables z\ , z 2 evolve between them in such a way 
that R{z\),R{z 2 ) remain non-negative. This corresponds to a quasiperiodic motion of 
the sphere. 

In the sequel we assume that the moments of inertia I± , 1% , I3 corresponds to a 
physical rigid body, i.e., that the triangular inequalities 7, + Ij > Ik are satisfies. For 
concreteness, assume also that d < Ji < J2 < J3. This also implies < A3 < A 2 < A\ 
and < ai < 02 < 03. 

As follows from the first expression in (|3.4[) . in the real case the factor G(zi, z%) is 
always positive. Hence, the right hand sides of (|3.2p are positive and the coordinates 
Hi are real and satisfy (n,n) = 1 if and only if z\ G [01,02] and z 2 G [02,03], like the 
usual sphcroconical coordinates. 

Next, we have 

Proposition 4.1. For the real motion, when the constants h, f are positive, and for 
any d, Ji satisfying the above inequalities, the roots c\ < C2 never coincide and 

2d — A 

ci < a± and if f jh — Ji, then c\ = ——. — t— < Qi, C2 = o-i- 

dAfrAj 

Proof. Set f /h = A G R. The roots Ci, C2 coincide with those of \E'(z)A + ip(z) and have 
the form 



Tr(AJ)(A + d)±VP 

Cl ' 2 = ddelA^ ' 

V = (Tr(AJ)) 2 (A + d) 2 - 8d det A(A + d) + 4d dct ATr J. 

The condition T> = gives a quadratic equation on A, whose determinant equals 
—d det I det A 2 , always a negative number. Hence V > and c\ < c 2 . 

Next, in view of (|4.ip . the condition c\ — a\ = also leads to a quadratic equation 
for A, again with always negative determinant. Then, evaluating c\ — a± for one value 
of A, we find c\ < a\ for any A 6 K. 

Finally, setting in (|4.1jl A = Ji and simplifying, we obtain the indicated above 
expressions for ci, C2. 



Combining the statement of Proposition 14. II with the permitted positions of z\,Z2, 
we conclude that, depending on value of C2, 

Zl G [ai,c 2 ], z 2 G [a 2 ,a 3 ], or Z\ G [ai,a 2 ], z 2 G [c 2 ,a 3 ]. (4.2) 

Then, in view of (|3.2|) , the vector n always fills a ring 1Z on the unit sphere S 2 — 
{(x,x) — 1} between the lines of its intersection with the cone 



Ji — d 



4— 1 Ji di- c 2 



= 0. 



Periodic solutions with bifurcations. As follows from Proposition 14. H the 
only periodic solutions with bifurcations can occur when the root c 2 coincides with 
ai,a 2 or aj^. This happens under the initial conditions 

uj t =ujj=0, n k = Q, (i,j, k) = (1, 2, 3), 



■2 



There is another type of periodic solutions corresponding to periodic windings of the 2-dimensional tori. 
However, the latter are not related to bifurcations and we do not consider them here. 
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when the sphere performs a periodic circular motion with = and one has (w, n) = 0, 
H = Jk^k, F\ — JfcU>%., which yields f /h = Jk- Then, in view of the above proposition, 
c 2 = Ofe and the polynomial R(z) in (|3.13j) has the double root Ofe, as expected. 

When A = f/h leaves the interval [Ji, J3], the root c 2 goes beyond of [ai, 0,3}. Then 
for R(zi), R(z2) to be both positive, one of Zi must violate the condition (|4.2[) . This 
implies that in the real case the quotient f/h belongs to [Ji, J3], and the bifurcation 
diagram on the plane (h, f) consists only of 3 rays f/h — J 1; J 2 , J3. 

Note that, according to the results of [5], a similar situation takes place for the 
Chaplygin sphere on the horizontal plane. 



The motion of the contact point on the fixed sphere. As mentioned 
above, in the generic case with F 2 = the contact point on the moving sphere given by 
the vector n belongs to the ring 71 on S 2 . 

Then the following natural question arises: does the contact point on the fixed sphere 
also belongs to a ring or it cover the whole sphere ? (Recall (|B|) that in the case of the 
Chaplygin sphere on a horizontal plane the contact point on the plane moves inside a 
strip, whose axis is orthogonal to the horizontal momentum vector.) 

To study the above problem we assume, without loss of generality, that the radius 
a of the fixed sphere is 1. Then the contact point on this sphere is given by the unit 
vector n as viewed in space. 

To describe the spatial evolution of n, introduce a fixed orthogonal frame O^rji^ with 
the center O in the center of the fixed sphere and the "vertical" axis 0( directed along 
the fixed momentum vector M. Then, in view of (|3.10p . the projection of n on 0( can 
be written in form 

" C = \W\ ^' M ^ = ^ aia2 ° 3 ( Jl _ J 2)(^2 - - Ji) n x n 2 n 3 

Z\ Z\ z 2 z 2 



(z\ - ai)0i - a 2 )(zi - a 3 )z 2 (z 2 - ai)(zi - a 2 ){zi - a 3 )zi 
which, following (|3.2[) and the expressions for zf,z 2 , after a simplification, reads 

7-J1J2J3 V (ci - zi)(ci - z 2 ) V( c 2 - zi)(c 2 - z 2 ) 



1 -1 

G(z\,z%) zi - z 2 



Z 2 Wi ZiW 2 



(zi - ci)(zi - c 2 ) {z 2 - ci){z 2 - c 2 ) 



(4.3) 



It follows that the right hand side of (|4.3[) is a quasiperiodic function of time. 
One can show that under the conditions (|4.2p and c\ < a\ < a 2 < a 3 , the function 



n^(zi,wi, z 2 ,w 2 ) is real and, regardless to signs of the roots w a — y / R(z a ), the function 
\uq \ has the absolute maximum in one of the vertices of the quadrangle Q C (zi, z 2 ) = M 2 



defined by (|4.2j) . In two other vertices of Q this function is zero. 

Calculating |n^| in the vertices of Q, we find that for c 2 ^ ai its maximum is strictly 
less than 1. It follows that the trajectory n(t) on the fixed sphere lies between the 
"horizontal" planes £ = v < 1. 

To describe the trajectory n(t) on the fixed sphere in the "longitudinal" direction, 
apart from the fixed momentum vector M it is good to know another fixed vector which 
can be expressed in terms of u>, n. However, it seems that such a vector does not exist, 
and for this reason we introduce the longitude angle ip between the axis 0£ and the 
vertical plane spanned by M and n. Introduce also the the longitudinal unit vector 
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u = M X n/\M x n\. Then wc find 



\M\ 



\M x n 



tt, — n 

dt 



\M\ 



(Mxn,jn) 
(M xn,M xn) 



where is the absolute derivative of n expressed in the coordinates of the moving 
frame. In view of the second vector equation in (2.1) with k = — 1, 



dt 



n = h + Lvxn = 2ojxn. 



Hence, we get 



Next, using the expressions 



4> 



\M\ 



(M x n, 2lj x n) 



2(uj x n)i = 
(M x n )i = 



n i (J l - d) 
G 

n i (J l - d) 
2G 



(M xn,M xn)' 
, (|3.10p , we obtain 

dA i A k z 2 + At - 2d . dA j A k zi+A i ~2d 



a, 1 z 1 



1 



-zi + 



a t 1 z 2 - 1 



(2JjJ k - dAi)z 2 - 1 . , (2J 3 J k - dAi)zi - 1 . 



K lz i - !) z 2 

= (1,2,3). 



(a,- ^2 - l)zi 



(4.4) 



Substituting these formulas into (|4.4[) and expressing the derivatives z a in terms of 
Zi, Wi, Z2, W2, one finds the derivative i/iasa symmetric function of (zi, w±) and (z2, W2), 
that is, as a quasiperiodic function of t. Its integration yields ip(t), which, together with 
(|4.3p . provides a complete description of the contact point on the fixed sphere. 



5 A special case of periodic motion. 

Apart from the particular case of the motion with F2 — 0, there is another special 
case, when this integral takes the maximal value, that is, when An is parallel to the 
momentum vector M. In this case M = hAn, h =const. In view of (|2.1jl . this implies 

Jo; — d(u>, n) = h An and co = hJ^ 1 ( An + —{An, J _1 n) J , (5.1) 

V p ) 

p being the same as in (|2.5[) . 

Substituting the expression for u> into the second equation in (|2.2p and simplifying, 
we get the following closed system for n: 

J\ + J2 + J?, — 2d x 
n = h — (n x J n) . (5.2) 

It has two independent integrals (n,n) and (n, J _1 n) or (An, An), which implies that 
the factor F is constant on the trajectories and that the system has the form of the Euler 
top equations. As a result, in the general case the components of n and w are expressed 
in terms of elliptic functions of the original time t and their evolution is periodic. 

This situation is similar to that of the special case of the motion of the Chaplygin 
sphere on a horizontal plane, when the momentum vector M is vertical, and when the 
solutions are elliptic in the original time. 
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6 Theta-function solutions in the case F2 = 0. 



In order to find explicit solutions for the components of u, M, n, and other variables, we 
first remind some necessary basic facts on the Jacobi inversion problem and its solution. 

Solving the Jacobi inversion problem by means of Wurzelfunktionen. 

Consider an odd- order genus g hyperelliptic Ricmann surface T obtained from the affinc 
curve 

= R(X)}, R(X) = (A - Ei) • • • (A - E 2g+1 )}, 

by adding one infinite point oo. Let us choose a canonical basis of cycles ai, . . . , O s , 
bi, . . . , b g on r such that 

^ o a, ; = bi o b 3 ■= 0, a,-b, <%, i,j = l,...,g, 

where 71 o 72 denotes the intersection index of the cycles 71 , 72 (For real branch points 
see an example in Figure RTTj) . Next, let u>i, . . . , a) g be the conjugated basis of normalized 
holomorphic differentials on T such that 

j> U>i = 27TJ 5ij , j = V^T- 

The g x g matrix of 6-periods Bij — § b uji is symmetric and has a negative definite real 
part. Consider the period lattice A = {27rjZ» + BZ 3 } of rank 2g in C 9 = (Z u . . . , Z g ). 
The complex torus Jac(r) = C 9 /A° is called the Jacobi variety (Jacobian) of the curve 
r. For a fixed point Po the Abel map 

^:r^Jac(r), A{P)= [ (wi,...,a) g ) T 

J Pa 

describes a natural embedding of the curve into its Jacobian. 




Figure 6.1: A canonical basis of cycles on the hyperelliptic curve represented as 2-fold covering of the 
complex plane A. The parts on the cycles on the lower sheet are shown by dashed lines. 

Now consider a generic divisor of points Pi — (Ai, fii), . . . , P g = (X g , on it, and 
the Abel-Jacobi mapping with a basepoint Po 

Pi fPo 

Q + ---+ lu = Z, (6.1) 

'Po J Pa 

(q u ...,lj b ) t , z = (Zi,...,z g f ec°. 
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Under the mapping, symmetric functions of the coordinates of the points Py, . . . ,P g are 
2g-fold periodic functions of the complex variables Z\ 1 . .. ,Z g with the above period 
lattice A (Abelian functions). 

Explicit expressions of such functions can be obtained by means of theta-functions 
on the universal covering C 9 = [Z\, . . . , Z g ) of the complex torus. Recall that customary 
Riemann's theta-function 8(Z\B) associated with the Riemann matrix B is defined by 
the serie^l 



9{Z\B)= J2 cx P ((BM,M) + (M,Z)), 



(6.2) 



MGZf 



(M, Z)=J2 M i Z i' ( BM > M ) = £ B ij M i M 3 



Equation 6(Z\B) = defines a codimension one subvariety S Jac(r) (for g > 2 
with singularities) called theta- divisor. 

We shall also use theta-functions with characteristics 



a= (ai, . . . , a g ), /3 = (/3i 



aj,/3j G 



which are obtained from 8(Z\B) by shifting the argument Z and multiplying by an 
exponen10: 



(Z) 



QL\ • ■ ■ OLg 



(Z) = exp{(Ba,a)/2 + (Z + 2irjf3,a)}6(Z + 2njl3 + Ba) 



All these functions enjoy the quadiperiodic property 



0. 



(Z + 2-kjK + BM) = exp(27rje) exp{-(BM, M)/2 - (M, Z)}9 



(Z), (6.3) 



(a,K)-{(3,M), 



Now for a generic divisor Pi — (Ai, /Zi), . . . , P g = {\ g .ji g ) on T, introduce the 
polynomial U(X,s) = (s — Ai) ■ • • (s — A g ), A e C. It is known (see e.g., [1] [2]) that 
given a generic constant C 7^ Ei, then under the Abel mapping ()6. 1(1 with Po = 00 the 
following relations hold 



tf(A,C) = (C-Ai)-..(C7-A fl ) = x 



g[A] (Z-g)g[A](Z + g) 
6*[A](Z) 



(6.4) 



q = A(C,VR(C)) = I (Qx,...,Q g ) T , 



where k is a constant depending on the periods of T only. 

These relation can be generalized in different ways as follows. 



3 The expression for 0(Z) we use here is different from that chosen in a series of books on theta-functions 
by multiplication of Z by a constant factor. 

4 Here and below we omit B in the theta-functional notation. 
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Theorem 6.1. (see, e.g., [TJ[2j[T6]). Under the Abel mapping (|6.ip with Pq — oo the 
following relations hold 

VuiKW) = s/iEi - Ai) ■ ■ ■ (Ei - X g ) = h 6[ ^^ ] , (6.5) 
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v Mfc y/U^EJy/U&Ej) 8{A + i hj }(Z) 

tiU ¥k (^-Xi) (Ek-XkXEj-Xk) 11 0[A](Z) ' [ -> 

i,j = l,---,2g- 1, i^j, 

where ki , kij are certain constants depending on the periods of T only, and 

are half-integer theta- characteristics such that 

27rj^'+B»jJ= / uj (mod A), (6.7) 



2mA" + BA' = K (mod A), and Vv=Vi+V 3 (mod Z 2ff ), 
/C G C 9 being the vector of the Riemann constants. 

Apparently, relations (|6.6[) were first obtained in the explicit form by Konigsberger 
(|22|). Earlier, expressions (|6.5p had been considered by K.Weierstrass as generalizations 
of the Jacobi elliptic functions sn(Z), cn(Z), and dn(Z). This set of remarkable relations 
between roots of certain functions on symmetric products of hyperelliptic curves and 
quotients of theta-functions with half-integer characteristics is historically referred to 
as Wurzelfunktionen (root functions). 

One can show (see, e.g., [131 H]) that for the chosen canonical basis of cycles 
<*i, . . • , a g , bi, .. . ,b g on T, 

A' = (1/2,...,1/2) T , A" = ( 5 /2,( 5 -l)/2,...,l,l/2) T (modi). (6.8) 

In the case g — 2, all the functions in ()6.5|) . (|6.6p are single- valued on the 16-fold 
covering T 2 -^Jac(r) with each of the four periods of Ao doubled, so that T 2 and Jac(r) 
are transformed to each other by the change Z — ► 2Z. In view of ()6.7|) . (|6.8p . one has 

(1/2 l/2\ A (0 l/2\ . / 1/2 



1/2 y ' 1 11 \0 1/2/' 1 u \l/2 1/2 

l'V2 0\ A/2 0\ A/2 1/2 N . 

" 3= (l/2 I/2J' A + ^ = (,l/2 0)' A + r ' 5= h/2 )' ,( '' ,! 



We shall also need the following modification of the Konigsberger formula (|6.6p . 
which, for our convenience, we adopt for the case g = 2. 

Theorem 6.2. Let C G C be a constant that does not coicide with Ei. Then under the 
Abel mapping (|6.ip with g — 2 and Pq —00: 



1 



Ai — A2 



Mi M2 



{E t - AiX^- - Ai)(C - Ai) (Ei - X 2 ){E 3 - A 2 )(C - A 2 ) 
f. 2 [A](Z) 9[A + mj ]{Z - g) 6[A + Vij ](Z + g) 
^ fl[A + ^ (Z) 9[A + W ] (Z) fl [A] (Z - ?) ff [A] (Z + g) ' 1 j 

i,j = !,•••, 2 5- 1, i^j, 
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where fey are constants depending on the periods of T only and, as in (|6.4[) , q = 
A(C, y/Wpj). 

Proof. Let us fix the point P2 — (A2,/i 2 ) in a generic position on the curve T and 
consider the following meromorphic function on this curve 

(A-_E,)(A-g 3 -)(A-C)) 
j,pj _ V + ^ 2 (A 2 -Si)(A2-i? 3 )(A2-C) 



(A-A 2 )(A-K i )(A-^)(A-C)- 



Due to the order of poles and zeros of z, w, and A — Ei, A — C on T, for any generic P2, 
the function /(P) has simple poles at P = P2, Pi, JSj, Q_, Q+ and does not have a pole 
neither at 1P2 = (A2, —^2), nor at any other point on I\ Next, /(P) has a double zero 
at co. Then, using the description of zeros of 0(Z), 9[A + r]i](Z), 6[A + r)ij](Z) one can 
show that up to a constant factor, /(P) = /(P) with 

2 [A](„4(P)-.4(P 2 )) 



'[A + mKAP) - AP2)) e[A + v 3 KAP) - a(p 2 )) 

g[A + gglQAQP) - g - A{P 2 )) 6{A + m] ](A(P) + g - A(P 2 )) 
e[A}(A(P) - q - A(P 2 )) e[A}(A(P) + q - A(P 2 )) 

Note that due to the quasiperiodic property of the theta-functions with characteristics, 
/(P) is a meromorphic function on G. 

Now setting P = iP\ = (Ai, — the function — /(P) transforms to the left hand 
side of (|6.10[) , and the argument A(P) — A(P 2 ) becomes —Z. Hence f(P) transforms 
to the right hand side of ()6.10|) . which proves the theorem. 



Combining Theorem 16.21 and formula (|6.4p . we obtain the following useful corollary. 
Proposition 6.3. Under the Abel mapping (|6.ip with g = 2 and Pq = 00, 

1 



Ai — A2 



(C-A 2 ) Ml (C-Ai)/* a 



(Ei - Ai)(P, - Ai) (Ei - X 2 )(Ej - A 2 ) 
6[A + r Hj ](Z-q) 0[A + 7fo](Z + q) 
f^A + ^K^A + ^KZ) 



= const, / 1 " Z;: 3l :Z u l\ lr7 71 1 ^ ■ ( 6 - U ) 



g = A(C, y/R(C)), i,j=l,...,2g-l, i^j, C £ E, t . 

Proof. Indeed, the left hand side of (|6.11|) is obtained from that of (I6.10|) by multipli- 
cation by (C — Ai)(C — A2), whose theta- function expression is given by formula (16. 4|) . 
Then the product of right hand sides of (|6.10p and (|6.4[) gives (|6.11|) . 

Explicit theta-function solutions. Now let T be the genus 2 curve 

{w 2 = R(z)}, R(z) = -(z - a 1 )(z - a 2 )(z - a 3 ) (z - c 1 )(z - c 2 ) , 

ci < C2 being the roots of f^f(z) + htjj(z). Thus we identify (without order) 

{Pi, • • • , P5} = {a.i,a 2l a 3 , c x ,c 2 }, 

and denote the corresponding half-integer characteristic r\i by r/ ai and r/ Ca . 

Next, choose the canonical basis of cycles as depicted in Fig. 16.11 and calculate the 
2x2 period matrix 

/ dz zdz 

Aij = (b ZUi, VJ\ — — 7==, H7 2 



y/Wz) ' ^R(z) 
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Then the normalized holomorphic differentials on T are 

E„ z J ~ 1 dz -, 

and the quadratures (|3.14p give 

"(zi.uji) r(z2,w 2 ) Azi.wi) p(z2,w 2 ) 

&1+ Wl = Zl, / (Il2 + / W2 = Z 2 , (6-12) 

^ oo ./ oo ./ oo 

Z 1 = 2C u t + Z 10 , Z 2 = 2C 21 t + Z 20 , (6.13) 

Zw,Z 2 o being constant phases. 

Now, comparing the last fraction in (13. 2|) with the expression (16. 5|) in Theorem 16.11 
wc find 



5^ ^-^-^ : = ^ + ^ , (m-,*:) = (1,2,3), (6.14) 
■ s /(a i -aj)(a i -ak) 9[/\\(Z) 

where Z = [Z\, Z 2 ) and the components of Z depend on r according to (|6.13p . 

To calculate fcj, we set here z\ — aj, z 2 = a^. Then, in view of (|6.12[) . the definition 
of 9 with characteristics and the quasiperiodic property (|6.3p . 

_ 9[A + 7 lai }(A(a J )+A(a k )) _ 0[A](O) fl[A + 77 ai ](0) 

0[A](„4(a,)+„4(a fe )) "SfA + ^KO)' ttotlS ' kl ~ 0[A](O) ' 

In order to express in theta- functions the factor G{z\, z 2 ) given by (|3.3|) . we fist note 
that it does not split into a product of linear functions in z\ and 22, hence one cannot 
use the formula (|6.4p . 

On the other hand, from the condition n\ + n 2 + n\ = 1 we find that 

_ dctl / Ji 2 J 2 q2 , ^3 q2 
G ~ dctJ ^ JT^ 1 J2 - d 2 + Ja^rf 1 
which, in view of (16.14[) . gives 



r- fdetT jnZ) ,,,,, 

VG(zi '^^VdetT?[Ap) ! ( } 

E(Z) = ^tf 2 [A + ^ 01 ](Z) + ^0 2 [A + r^KZ) + ^tf 2 [A + Vas ](Z). 

^1 J2 ^3 

As a local singularity analysis shows, in general the function T,(Z) has zeros of first 
order, hence it cannot be a full square of another theta- function expression. 

Now, using theta- function expressions (16.14|l . (16.15P in formulas (13. 2p . (j3 . 1 5[> . (|3.16p 
and applying also the Wurzelfunktionen (|6.6p . (|6.1ip with C = 1/(gL4,), we arrive at 
the following theorem. 

Theorem 6.4. The generic theta-function solutions for the Chaplygin sphere-sphere 
problem in the case F 2 = have the form 

ni(T)=Ki vm ' ( } 

M t (r)= Vi 9[A + r ^" m , (6.17) 

f , OJA + V aj + Va k }(Z - gi) g[A + ^ + ria k }(Z + gQ 
usAt) = £i , , 6.18 

e[A](z) ■ jnzi 

q. t = A(l/(dAi),y/R(l/(dAi))), Ki,Vi,ei = const, (i,j,k) = (1,2,3), 
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where the characteristics are given in (|6.9p and Z\ , Z2 depend linearly on r as described 

m mm). 

We do not give explicit expressions for the constants Ki,Vi,£i here. 

Note that due to presence of the square root, the variables rij, Mi,uii are not mero- 
morphic functions of Z\,Z2 and therefore, of the new time r. This stays in contrast 
with the solutions of the classical Chaplygin sphere problem, which, after a similar time 
reparameterization, become meromorphic (see [111 117) ). 

Next, comparing the expression (|3.17|) with the Wurzelfunktion (|6.6|) and (|4.3|) with 
(|6.1ip . assuming C = 0, we also find 

9[A + Vci +n C2 ](Z) 

M = v — 9[Am — ' ( } 

g[A + r /ci + n C2 ]{Z - q) 0[A + n Cl + n C2 }(Z + g) 
n C — 6 ; ? (o.zUI 

6[A](z) vnz) 

q = A(0, \J R{0)), v,Q = const. 

The second formula, together with (|6 . 13|) . describes the altitude of the contact point on 
the fixed sphere as a function of r. 

Finally, given the expression (|6.15[) for the factor G, the original time t can be found 
as a function of r by integrating the quadrature J3.12[) . 

Appendix. Separation of variables via reduction to a 
Hamiltonian system 

As mentioned above, the substitution (|3.2p is quite non-trivial and can hardly be guessed 
a priori. Below we describe how one can obtain it in a systematic way. 

Al. Reduction to a Hamiltonian system on S 2 

First introduce the spheroconical coordinates u, v on the Poisson sphere {n,n) = 1: 

"•'w^m^k- '*'* k *'- J, = I,+i (6 - 21) 

Then, under the substitution (|3.1[) . the equations (|2.2p with k = — 1 give rise to the 
following Chaplygin-type system on S 2 

ddT dT , x ddT dT 
at ou ou at ov ov 

T = ^{b uu u 2 + b uv iiv + b vv v 2 ), $ = (a u ii + a v v), 



where T is the energy integral (|2.3|) expressed in the spheroconical coordinates under 
the condition F2 = and $ is linear homogeneous in it, v. Explicit expressions for the 
coefficients of T, $ are quite tedious, so we do not give them here. 

dT dT 

Introducing the momenta P u = — — , P v — — — , this system can be transformed to a 

ou di> 

Hamiltonian form with extra terms, which possesses invariant measure N du dv dP u dP v 
with the density 

2uv +(u + v)(2d + o-i) + a 2 ~ rigi , 3 2 

A* = ^ (4a 3 + 2ct\ct2 — Q-i — da 1 + 

Ay/det(J — dn®n) 

+(af - 2a 2 + 4dctx)(u + v) - 4d(u + v) 2 )~\ 
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where a\ = J2 Jh a 2 — J2 J?i a 3 = Ji-h-h- 

According to the Chaplygin theory of reducing multiplier ([8]), after the time repa- 
rameterization N(u, v) dt — dr the system (|6.22|) is transformed to the Lagrange form 



d dT dT _ d dT dT _ 
dr du' du rl ~ r ^'"' 



du 



dv 



dr dv' dv 



dr ' dr 



(6.23) 



As a result, under the time reparameterization we obtain an intcgrable Hamiltonian 
system on the cotangent bundle T*S 2 with local coordinates u,v,p u — dT / du' ,p v = 
dT/dv'. 



A2. Separation of variables 



Equations (|6.23|) possess 2 homogeneous quadratic integrals, which come from H, Fi 
(12.31) and which can be written in the form 



H = T = ^{g uu {u'f + 2g uv u'v' + g vv (v') 2 ), 

Fx = \{G uu {u'? + 2G uv u'v' + G vv {v') 2 ). 

Explicit expressions for the coefficients g uu , . . . , G vv as functions of u, v are suppressed 
due to their complexity. 

According to the result of Eisenhart [12] (see its modern accounting in [24]). sepa- 
rating variables S\,S2 can be chosen as the roots of the equation 



det(G - sg) = 0, 



(6.24) 





G U v 




Quu Quv 




Gyy 


,g = 


Quv Qvv 



G = 



Note that the roots depend only on the local coordinates u, v on the configuration 
space S 2 . 

The spheroconical coordinates (|6.2ip depend explicitly on the roots s±, S2 as follows 



1 



1 



where 



±s 2 Q{si) ± s\Q(s 2 ) 
4d(si - s 2 ) 



y 



±Q( Sl )±Q(s 2 ) 



2d(si 



S2 



Q(s) = \lb lS 2 -b 2 s + bl 

2 



— (Tr(JA)) 2 - -ddet A 
16 y V " 2 



2detJ- -dTr(JA), 



b 2 = det JTr(JA) - id(Tr(JA)) 2 - ^ddet ATrJ + d 2 det A. 

dT 

In the new variables si,s 2 and the conjugated momenta p\ = -tttiPz 
integrals take the Liouville form 



(6.25) 
(6.26) 
(6.27) 
dT 



ds' 2 



the 



Si - s 2 



82{S2) -Vl F\ 



si - s 2 



fgjjXfl) 2 
Pi 

si - s 2 



SlS 2 (s 2 ) 2 
Si - s 2 



(6.28) 
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where 



_ 2(8a: 3 + 8{d - e)x 2 + (2e 2 /3 - Ade)x - 4 7 - d/3 + y^Mx)) 
X ' ~ 7 (2x - e + 2d) 2 ' 

A = x 2 {0 2 + 8ad) + 2a;(4/?7 + d/3 2 - 2dae + 4ad 2 ) + (4 7 + d/3) 2 , 



and we used the notation 

a=(J 2 + J 1 - J 3 )(-J 2 + J2- + -h + h) = -A^Az, 

P = Jf + J 2 + J| - 2Ji J 2 - 2J 2 J 3 - 2J 3 Ji = -A XJ 4 2 - A 2 A 3 - A 3 A U (6.30) 

7 = JiJ 2 J 3 , e = Ji+J 2 + J 3 . 

Due to the Hamilton equations with the Hamiltonian H in (I6.28|) . the evolution of 
si,s 2 is descibed as follows 

„/ _ VVi yi j _ \fih V2 , g 31 , 



2si — e + 2d si — s 2 2s 2 — e + 2d s 2 — s\ 

Vi = y/Y(xi), 



Y(x) = A(x) ■ (hx - f) 8x 3 + 8(d - e)x 2 + (2e 2 - /3 - 4de)x - c + y/A(x) , (6.32) 

where h, f are the constants of the integrals H, F\. 

Hence, we performed a separation of variables, however the evolution equations 
(|6.3ip have a quite tedious form. 

One can show that the equation y 2 = Y(x) defines an algebraic curve of genus 2 
on the plane C 2 = (x,y). According to the theory of algebraic curves, any curve of 
genus 2 is hyperelliptic and can be transformed to a canonical Weierstrass form by an 
appropriate birational transformation of the coordinates x, y. 

One of such transformations is induced by the chain of substitutions x — > £ — > z 

46 3 £ -(4detJ-dTr(AJ))z + TrJ-2d 



{i + b 2 ) 2 -4&i& 2 ' s 2detAdetI 

h,b 2 ,b 3 being defined in (j6~26"]) . (f6T27|) . It converts A (a;) in ([6T32]) . as well as Q(x) (j6?25|) 
into full squares. 

After some tedious calculations, one finds that in the new variables z\ = z(x\), z 2 = 
z(x 2 ) the expressions (|6.21l) take the form (|3.2|) . which ensures the reduction to hyper- 
elliptic quadratures in the canonical form (|3.14|) . 
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